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基于 二 次 插值 的 天 牛 须 搜索 算法 
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(江西 理工 大 学 , 信息 工程 学 院 , 江西 赣州 341000) 


摘 Xe: 针对 天 牛 须 搜索 算法 在 高 维 空间 中 搜索 精度 低 和 易 陷 入 局 部 最 优 的 问题 进行 了 研究 ， 提 出 一 种 新 的 天 牛 须 
优化 算法 -基于 二 次 插值 的 天 牛 须 搜索 算法 , 称 之 为 QIBAS。 算 法 在 天 牛 进行 移动 后 , 将 天 和 牛 当 前 位 置 左右 两 触须 作 
为 插值 坐标 点 ， 利 用 二 次 插值 生成 一 个 新 的 解 ， 再 对 比 插 值 产生 的 解 与 当前 最 优 解 、 全 局 最 优 解 的 适应 度 值 ， 更 新 
全 局 最 优 解 。 对 多 个 单 峰 函 数 和 多 峰 函 数 进行 数值 仿真 测试 ， 其 维度 分 别 取 100、500、1000、5000、10000， 仿 真 结 
果 表 明 ， 引 入 二 次 桂 值 有 效 提升 了 7 BAS 算法 跳出 局 部 最 优 的 能 力 。QIBAS 在 求解 最 优 值 时 ， 其 求解 精度 有 极 大 的 
提升 ， 收 么 速度 也 有 较 明 显 提升 ， 改 进 算 法 的 有 效 性 得 以 验证 。 
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Beetle antennae search based on quadratic interpolation 


Liao Liefai, Ouyang Zongying 
(Jiangxi University of Science & Technology, School of Information Engineering, Ganzhou Jiangxi 341000, China) 


Abstract: Focus on the problems that low-precision and easily to fall into local optimum result in the high-dimensional space 
by the longhorn search algorithm, we proposed a new optimization algorithm of the short-lived beetles-the short-lived search 
algorithm based on quadratic interpolation which is called QIBAS. The algorithm takes tentacle's left and right tentacles as 
interpolation coordinate points after the beetle moved. And generate new solution with a second interpolation, then compare 
to the current optimal solution and the global optimal solution with the fitness of interpolation solution. At the same time to 
update the global optimal solution. Numerical simulation tests were performed on multiple unimodal and multimodal functions, 
and their dimensions take as 100, 500, 1000, 5000 and 10000. The simulation results show that the introduction of quadratic 
interpolation effectively improves the ability of the BAS algorithm to jump out of the local optimum. When QIBAS solves 
the optimal value, its solution accuracy is greatly improved, and the convergence speed is also significantly improved. The 
effectiveness of the improved algorithm is verified. 
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了 BAS 随机 种 子 导 致 的 不 同 优化 结果 和 优化 精度 低 的 问题 。 
文献 [17] 将 模拟 退火 的 蒙特 卡 洛 准则 引入 BAS 中 ， 加 快 算法 
能 优化 算法 是 人 类 对 大 自然 的 探索 而 形成 的 ， 对 自 跳出 局 部 最 优 ， 并 将 优化 后 的 算法 应 用 于 分 布 式 电源 选 址 定 
然 界 各 类 生物 行为 的 仿真 模拟 。 智 能 算法 广泛 应 用 在 图 像 处 容 ; 文献 [18] 将 BAS 与 粒子 群 算法 相 结合 ， 提 升 了 其 收敛 速 
理 、 工 业 控 制 、 生 产 调 度 、 模 式 识 别 等 诸多 领域 。 智 能 算法 度 和 搜索 精度 。 文 献 [19] 结 合 BAS 和 梯度 下 降 ， 选 择 梯 度 下 
的 种 类 也 颇 多 ， 有 蚁 群 算法 中 、 粒 子 群 算法 中、 遗传 算法 中、 降 方向 或 天 牛 “ 左 右 须 ”最 优 解 方向 进行 迭代 。 
人 工蜂 群 算法 内、 鲸鱼 算法 本 等 等 。 在 处 理 优 化 问题 时 ， 这 些 上 述 对 BAS 算法 的 改进 均 取 得 一 定 效果 ， 但 对 于 BAS 
算法 都 有 一 定 缺 陷 , 如 蚁 群 算法 所 需 设置 的 参数 较为 复杂 , 且 算法 在 高 维 空间 中 搜索 精度 低 的 情况 ， 上 述 算法 均 未 提出 解决 
计算 量 大 ;粒子 群 算法 在 处 理 问题 时 易 早 熟 收 全 ， 尤 其 是 复 方案 。 经 研究 发 现 ， 在 高 维 空间 中 甚至 出 现 天 牛 个 体 不 移动 的 
杂 多 峰 问 题 ， 遗 传 算法 需 对 问题 进行 编码 和 解码 ， 其 复杂 度 情况 ， 这 大 大 降低 了 算法 的 寻 优 精度 。 而 针对 该 问题 ， 本 文 将 
较 高 。 国内 外 众多 学 者 也 针对 这 些 问 题 对 算法 进行 优化 [10。 二 次 插值 引入 到 BAS 算法 中 。 在 天 牛 进行 一 次 移动 后 ， 借 助 天 
天 牛 须 搜索 算法 (beetle antennae search, BAS)! Jiang 牛 的 左右 两 须 ， 利 用 二 次 插值 法 产生 新 的 解 ， 再 对 比 新 解 与 当 
和 Li Æ 2017 年 在 提出 的 新 型 智能 优化 算法 。 该 算法 模拟 天 前 全 局 最 优 解 ， 以 此 更 新 全 局 最 优 解 。 仿 真 结果 表明 ， 改 进 的 
牛 竟 食 行 为 ， 以 此 实现 对 复杂 优化 问题 的 寻 优 求解 。BAS 算 算法 (QIBAS) 在 求解 高 维 问题 时 相 较 于 天 牛 须 算法 (BAS)、 粒 子 
法 广泛 应 用 于 各 领域 。 如 文献 [13] 将 BAS 算法 应 用 于 微 电 网 群 算法 (PSO) 和 遗传 算法 (GA), 在 求解 精度 上 有 了 进一步 的 提升 。 
多 目标 能 量 管理 中 ， 为 微 电 网 的 能 源 优 化 管理 和 调度 提供 了 TN 
新 的 解决 方案 ， 文 献 [14] 用 BAS 优化 神经 网 络 权 值 和 阔 值 ， 1 ”天 和 牛 须 搜索 算法 
来 预测 喷射 灌浆 混凝土 的 无 侧 限 抗 压强 度 。 但 针对 BAS 算法 天 和 牛 须 搜索 算法 区 别 于 其 他 优化 算法 ， 该 算法 只 包 
易 陷 入 局 部 最 优 、 稳 定性 不 足以 及 仅 能 解决 单 目标 优化 问题 个 个 体 ， 不 需要 清楚 具体 函数 和 梯度 信息 就 能 进行 寻 优 ， 基 
5 的 缺点 ， 学 者 们 提出 了 一 系列 的 改进 方法 。 如 文献 [16] 受 ” 此 该 算法 运算 量 小 、 寻 优 速 度 快 且 易 于 实现 。 
群体 优化 算法 启发 , 将 群体 引入 BAS 算法 中 , 提出 了 天 和 牛 群 1.1 基本 原理 
优化 算法 (beetle swarm antennae search, BSAS), BSAS 克服 天 牛 须 搜索 算法 的 仿生 机 理 是 ， 天 牛 在 竟 食 时 不 清楚 食 
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的 过 程 。 


天 牛 须 搜索 算法 的 智能 机 理 是 ， 


右 移动 ， 最 后 寻找 到 食物 的 确切 位 
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图 1 天 牛 搜索 食物 过 程 
Fig. 1 Beetle searching for food 


。 图 1 为 天 牛 搜索 食物 
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单个 天 牛 在 搜索 食物 时 ， 


食物 气味 浓度 即 为 适应 度 函 数 ， 这 个 函数 在 空 


有 不 同 的 值 ， 
对 比 二 者 气味 浓度 ， 


m 


mAH 


探测 的 气味 浓度 相差 
DAIEN 


如 下 几 个 步骤 : 
a) 天 牛 初始 朝向 创建 


其 中 ranas) 为 随机 函数 ， 为 空 


_ rands(K,1) 
|| rands(K,1) || 


间 维 度 。 


al 


F 须 算法 的 寻 优 过 程 。 


间 的 不 同位 置 


两 须 分 别 采 集 自身 附近 两 点 的 气味 值 ， 

取 最 优 值 ， 如 此 过 代 ， 知 道 左 右 两 须 所 
是 一 定 精度 则 停止 兴 代 ， 所 找到 的 中 
司 食物 浓度 最 高 点 ， 即 食物 所 在 地 。 

天 牛 搜寻 食物 的 过 程 就 是 天 4 


分 为 


(1) 


录用 定稿 BINE, F: 基于 二 次 插值 的 天 牛 须 搜索 算法 第 38 卷 第 3 期 
物 的 具体 位 置 ， 利 用 头 上 一 对 触须 来 感应 食物 气味 的 浓度 ， 有 跳出 局 部 最 优 的 能 力 , 因此 BAS 算法 在 每 次 寻 优 的 结果 都 
如 果 左 触须 感受 到 的 气味 浓度 高 ， 天 牛 朝 左 移动 ， 反 之 则 朝 有 不 同 。 针 对 高 维 问题 ， 由 于 BAS 算法 是 单个 体 搜索 算法 ， 


单个 天 牛 对 食物 气味 浓度 的 分 辩 能 力 不 足 以 让 其 在 高 维 空间 
中 找到 最 优 解 ， 这 就 可 能 导致 天 牛 位 置 不 更 新 ， 完 全 陷入 
部 最 优 ， 大 大 降低 了 BAS 在 高 维 空间 的 搜索 精度 。 
2 天 和 牛 须 搜索 算法 的 改进 
2.1 算法 改进 

本 文 为 解决 BAS 在 高 维 空间 中 优化 精度 低 的 情况 , 引入 
了 二 次 插值 算 子 式 (6)。 Wd ds 种 用 于 在 确定 初始 区 
间 中 搜索 极 值 点 的 方法 ， 种 曲线 拟 合 方法 。 其 利用 目标 
函数 在 若干 点 的 信 ， 息 构成 一 个 与 目标 函数 值 相近 的 低 次 多 项 
式 ， 再 用 此 式 的 最 优 解 作为 函数 的 近似 最 优 解 ， 随 着 区 间 的 
逐步 缩短 ， 多 项 式 的 最 优 解 与 原 函 数 的 最 优点 之 间 的 距离 逐 
步 减 小 ,最 后 满足 一 定 精度 要 求 为 止 。 其 原理 图 如 图 2 所 示 。 


a 


—— 


图 2 
Fig.2 Principle of quadratic interpolation 


二 次 插值 原理 


本 文 在 每 一 次 天 牛 位 置 更 新 时 ， 将 左右 两 须 所 在 的 点 作 
为 二 次 插值 的 区 间 ， 构 成 一 个 与 目标 函数 值 相近 的 低 次 多 项 
式 ， 利 用 该 式 求 出 最 优 解 ， 根 据 最 优 解 再 缩小 区 间 ， 迁 代 进 
行 ， 最 后 再 满足 一 定 精度 的 情况 下 ， 得 到 一 个 区 间 最 优 解 ， 
间 最 优 位 置 。 并 与 全 局 最 优 解 进行 对 比 ， 取 两 者 最 优 为 
新 全 局 最 优 解 。 利 用 这 种 方法 大 大 提升 了 天 牛 须 搜索 算法 的 
局 部 搜索 能 力 ， 进 而 提高 了 天 牛 在 高 维 空间 中 跳出 局 部 最 优 


的 能 力 、 搜 索 精 度 和 收敛 性 能 

设 天 牛 的 第 1+ 代 在 高 维 空间 中 质心 位 置 为 
X,2Quxoesxx ， 当 前 天 牛 两 须 位 置 坐标 分 别 为 
X, = nX eA) REX, =x,%2,…Xx) ， 当 前 全 局 最 优 位 置 为 


X, = (Ki, poss 


wx)”， 则 二 次 插值 生成 的 位 置 为 

l — C ox) f G5) OG. 7 X3) f 08) t GG — Xi) f O9) 

2 [Co — 25) f Q5) + o — x4 f 09) + Qo — X) f 09,)] + eps 
kzLl2,..K (6) 


Xi 


设置 步 长 因子 5 。 天 牛 的 区 域 搜索 能 力 由 < 决定 ， 为 了 
提高 天 牛 的 搜索 范围 ， 本 文选 取 了 较 大 初始 步 长 ， 并 线性 方 
式 递减 步 长 。 公 式 如 下 ; 

&a-7ó*eta (t=1,.…,n) (2) 

其 中 eta 为 递减 因子 ， 范 围 为 [0.] ，1 为 当前 迭代 次 数 ，n 为 

总 迭代 次 数 。 

b) 建立 天 和 牛 两 须 位 置 坐标 
— f 

2 Qs) 6) 
Tetti 

其 中 : 尺 表 示 左 须 的 位 置 坐 标 ; 表示 左 须 的 位 置 坐标 ; x 

表示 人 迭代 次 数 为 ! 时 天 牛 的 质心 坐标 ;ww 为 两 触须 之 间 长 度 ， 

其 值 应 足够 大 以 覆盖 适当 的 搜索 区 域 ， 以 便 在 开始 时 跳出 局 


部 最 小 点 。 


计算 适应 度 函 数 f0)， 


FO) 和 f(x,) 


fO 指 的 是 左右 须 气味 浓度 ， 即 


的 大 小 。 


的 选取 ， 将 在 后 文 详细 阐述 。 


fitness =T X (d, da) 


£O 其 根据 实际 需求 定 ， 对 


适应 度 函 数 


(4) 


Hp: 70 为 适应 度 函 数 ;，eps 是 个 非常 小 的 正 数 ， 为 防止 分 
母 为 0。 
QIBAS 算法 的 具体 步骤 如 下 : 


a) 初始 化 天 牛 。 根 据 式 (1) 随 机 生成 天 牛 的 初始 朝向 


d,-(d.d,--d,Y ， 随 机 生成 天 牛 的 初始 质心 位 置 
X, =ar X), 并 根据 式 (4) 求 得 适应 度 值 为 £O ,全 局 最 


优 值 Xi = Xo , 适应 度 值 为 1(Xs)， 设 定 步 长 因子 为 ,递减 


其 中 : da 为 第 i 个 样本 模型 输出 值 ，4, 为 第 i 个 样本 实际 值 。 


更 新 天 牛 的 


fo)» f(x), 
如 下 : 


其 中 表示 1 次 迭代 步 长 因子 ; 
性 能 分 析 


1.2 
天 牛 须 ] 


立 置 ， 比 较 左右 触须 适应 度 值 的 大 小 ， 若 
天 牛 向 左 移动 ， 反 之 ， 向 右 移动 。 其 更 新 公式 


x =x é *d* sign(f (x, —x,)) 


(5) 


sign() 为 符号 函数 。 


搜索 算法 的 优化 过 程 
别 


左右 两 端的 食 
浓郁 则 天 牛 就 向 哪 边 移动 ， 
提升 寻 优 速度 ， 但 也 会 导致 天 牛 陷入 局 部 最 优 ，f 


VI RET H 


3E 


时 ， 主 要 是 依靠 天 牛 
， 哪 边 触须 接收 到 的 食物 气味 
虽 能 大 大 加 快 天 牛 位 置 的 


触须 对 


P 


更 新 和 


天 牛 又 没 


HF ea, 4t-0., 
b) 计算 天 和 牛 两 须 坐 标 和 适应 度 值 。 根 据 当 前 天 牛 的 位 置 
确定 天 牛 两 须 的 坐标 ， 如 式 (3) 所 示 ， 计 算出 左右 触须 位 置 坐 
F2 BEA X! = Oha) RI Xr m Ones xx t=1,2 n, K 
中 并 为 总 迭代 次 数 。 根 据 式 (4) 求 得 两 须 的 适应 度 值 £OCO 和 
fO. 

c) 二 次 插值 和 天 牛 位 置 更 新 。 根 据 式 (5) 来 天 牛 移动 的 位 
Ti^ X, ,并 求 出 此 时 质心 位 置 对 应 的 适应 度 值 f(X,) ,根据 
式 (6), 求 得 插值 产生 的 位 置 坐标 为 X? = Gh" , EVER 
出 适应 度 值 为 7(X?) o HE SADO MIA), ESX) 则 
新 天 牛 位 置 坐标 ， 反 之 ， 则 保持 当前 位 置 


Hg E 
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d) 更 新 全 局 最 优 值 Xi 。 对 比 SAD 和 f(Xss) ， 取 两 者 算法 性 能 最 优 。 最 后 本 文 取 eta =0.1。 
中 最 优 值 替 换 全 局 最 优 值 ， 更 新 /Ce ) 和 Xow。 表 4 递减 因子 的 对 算法 性 能 影响 
e) 若 达 到 算法 终止 条 件 ， 则 停止 ， 否 则 +t=t+1， 返 回 步 Tab.4 Effect of decreasing factor on algorithm performance 
BE b). 函数 递减 因子 “最 优 解 “” 最 差 解 平均 收敛 代数 
2.2 ”改进 算法 时 间 复 杂 度 分 析 0.1 0 0 11.4 
本 文 的 改进 主要 是 在 天 牛 位 置 更 新 后 ， 利 用 二 次 插值 产 0.2 0 0 13.8 
生 一 个 异 于 当前 位 置 的 点 ,最 后 对 比 当前 位 置 , 取 两 者 优 值 。 0.3 0 0 16.5 
优化 问题 维度 为 XK 时， 本 文 算法 计算 复杂 度 分 析 如 下 : 初始 04 0 0 20.5 
化 天 牛 的 计算 复杂 度 为 0(K) ; 迭代 过 程 ， 计 算 两 须 坐 标 和 适 A 0.5 0 0 243 
应 度 值 ， 计 算 复 杂 度 为 0QK) ， 更 新 天 牛 位 置 和 计算 适应 度 EN 0.6 0 0 311 
值 , 计算 复杂 度 为 0QK) ,二 次 插值 产生 新 位 置 , 计算 复杂 度 0.7 0 0 42.6 
为 CCK) 。 每 一 次 迭代 的 计算 复杂 度 为 C6K) , 高 于 标准 天 牛 0.8 0 0 59.8 
须 算 法 的 OCK) 。 计 算 复 杂 度 可 近似 为 OUO 。 表 1 列举 了 遗 0.9 0 0 116.4 
传 算法 (GA)、 粒 子 群 算法 (PSO)、 标 准 天 牛 须 算法 (BAS) 及 本 1.0 0 0 807.5 
文 改进 算法 (QIBAS) 的 计算 复杂 度 。 表 中 为 种 群 数 0.1 10 10 8.2 
Al 算法 复杂 度 0.2 10 10 10.4 
Tab. 1 Algorithm complexity 0.3 9.99 10 12.5 
算法 计算 复杂 度 0.4 9.95 9.99 16.2 
GA O(NK + N?) 0.5 9.62 9.98 20.3 
Rosenbrock 
PSO O(NK +N) 0.6 9.92 10.0 26.1 
BAS O(K) 0.7 7.68 9.88 36.5 
QIBAS O(K) 0.8 5.12 9.62 57.7 
从 上 表 可 以 看 出 ， 算 法 计算 复杂 度 与 问题 维度 和 种 群 数 0.9 4.04 8.23 130.6 
量 有 关 ， 根据 大 0 的 概念 ， 若 优化 问题 维度 K 远大 于 种 群 数 1.0 4.97 7.97 853.2 
HN, GA 和 PSO 两 者 的 计算 复杂 度 均 为 OWK)， 而 BAS 0.1 -90 -90 8.5 
和 QIBAS 的 计算 复杂 度 为 0(K) 。 0.2 -90 -90 11.3 
2.3 ”改进 算法 参数 分 析 0.3 -90 -90 14.5 
根据 上 述 算法 模型 可 知 ， 改 进 天 牛 须 搜索 算法 需要 设 定 0.4 -90 -90 17.6 
的 参数 有 : 初始 步 长 印 ， 步 长 递减 因子 eta, WARKA N , e 0.5 -90 -90 21.7 
空间 维度 K 。 对 于 绝 大 多 数 应 用 问题 可 根据 需要 设 定 初始 步 0.6 -90 -90 28.9 
长 各， 本 文 取 值 为 1。 通 过 单 因 素数 值 实验 测试 函数 分 析 改 0.7 -90 -90 38.4 
进 算法 其 他 3 个 参数 对 算法 性 能 的 影响 ,测试 函数 如 表 2 所 示 。 0.8 -90 -90 59.8 
表 2 测试 函数 0.9 -90 -90 108.6 
Tab.2 Test function 1.0 -86 -81 337.6 
函数 名 称 表达 式 b) 最 大 迭代 次 数 对 算法 性 能 的 影响 
Sphere 2x TICK NUUCBLAT HHX n —1000,2000,3000,4000 。 对 测试 函数 
Rosenbrock 311006, ^x «Q7 x2] 进行 独立 20 次 数值 实验 。 实 验 结果 如 表 5 所 示 。 
Rastrigin ”> pe -10cos 27x, +10] 表 5 最 大 迭代 次 数 的 对 算法 性 能 影响 
算法 参数 的 经 验 设置 DE 3， 在 进行 单 因 素 实验 时 ， 固 Tab.5 Effect of maximum iterations on algorithm performance 
定 其 他 参数 ， 改 变 一 个 参数 进行 仿真 实验 。 函数 迭代 次 数 ” 最 优 解 ”最 差 解 。 ”平均 收敛 代数 
表 3 经 验 参数 设置 1000 0 0 116.2 
Tab.3 Experience parameter setting Sphere 2000 0 0 115.5 
算法 参数 变量 名 取 值 3000 0 0 117.6 
步 长 递减 因子 eta 0.9 4000 0 0 116.9 
迭代 次 数 n 1000 1000 4.96 8.01 30.6 
空间 维度 K 10 2000 4.95 7.89 30.8 
Rosenbrock 
a) 步 长 递减 因子 对 算法 性 能 的 影响 3000 5.28 8.12 30.5 
根据 常用 递减 因子 取 值 区 间 ee=[00 ， 设 置 步 长 因子 为 4000 5.06 7.88 31.5 
eta=0.1,0.2,0.3,0.4,…0.9,1。 对 测试 函数 进行 独立 20 次 数值 实验 。 1000 -90 -90 29.9 
实验 结果 如 表 4 所 示 。 人 2000 -90 -90 30.1 
步 长 递减 因子 的 设 定 是 为 了 算法 在 不 同时 期 有 不 同 的 搜 3000 -90 -90 29.6 
索 半 径 ， 在 算法 前 期 搜索 范围 很 大 ， 需 要 更 大 的 搜索 半径 ， 4000 -90 -90 29.9 
用 来 加 快 收敛 速度 。 随 着 迭代 次 数 的 不 断 增 加 ， 求 解 值 也 不 由 上 表 可 知 ， 最 大 迭代 次 数 对 于 本 文 算 法 的 影响 不 大 ， 
断 有 逼近 最 优 值 ， 此 时 搜索 半径 不 需要 太 大 ， 需 要 更 小 的 搜索 综合 考虑 本 文 取 n 22000 。 
范围 来 确定 最 优 值 。 — 
而 根据 表 4 得 出 ,对 Sphere 函数 而 言 当 递减 因子 ea=0.1 3 ， 仿真 测试 与 结果 分 析 
时 算法 性 能 最 优 . 对 了 osenbrock 函数 而 言 当 递减 因子 eta — 0.9 31 参数 设置 
时 算法 性 能 最 优 . 对 Rastrigin 函数 而 言 当 递减 因子 eta — 0.1 BT 本 文 基于 Matlab2018b 平台 仿真 测试 并 分 析 QIBAS 算 
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法 对 测试 函数 的 计算 性 能 。BAS 和 QIBAS 的 初始 步 长 乌 设 ”还 是 较 差 的 , 但 比 GA 算法 效果 更 优 。 在 与 标准 BAS 算法 比 
为 1， 递减 因子 分 别 取 为 eta, =09, eta, =01， 总 迭代 次 数 为 较 时 ， 除 了 太 不 如 三 个 函数 的 计算 结果 相差 不 大 ， 其 余 函 
n=2000 , 优化 问题 维度 分 别 设 为 Kk=100, 500, 1000, 共 运 行 50 ”” 数 计算 结果 的 平均 值 和 标准 差 都 较 优 。 由 此 可 以 看 出 改进 天 
次 测试 。GA 和 PSO 的 种 群 数量 均 设 为 200，PSO 中 学 习 因 ”人 牛 须 搜索 算法 的 计算 精度 和 稳定 性 都 有 较 明 显 的 提升 ， 验 证 


Rs 


子 均 设 为 1。 了 改进 算法 寻 优 的 有 效 性 。 
3.2 测试 函数 KO 单 峰 函数 
为 验证 算法 的 性 能 , 本 文选 取 单 峰 和 多 峰 等 共 11 个 标准 Tab.6 Unimodal function 
测试 函数 P9 进 行 测试 ， 测 试 函数 如 表 6. 7. 函数 测试 函数 函数 测试 函数 
3.3 实验 结果 与 分 析 fi Sphere Js Rosenbrock 
将 QIBAS 算法 与 GA ik. PSO 算法 以 及 BAS 算法 进 Schwefel Problem 2.22 fo Step 
行 对 比 测试 ， 各 运行 50 次 ， 取 平均 值 ， 并 计算 标准 差 ， 实验 f Schwefel Problem 1.2 f£ Quartic 
结果 如 表 8( 维 度 K 2200). 同时 ， xHT S 结果 进行 Friedman fa Schwefel Problem 2.21 fs Schwefel Problem 2.26 
非 参数 统计 分 析 ， 其 中 显著 性 水 平 设 为 5% ， 实 验 结果 如 表 9 表 7 多 峰 函 数 
所 示 。 表 10 是 各 算法 运行 50 次 太刀 三 个 函数 所 需 的 时 间 。 Tab.7 Multimodal function 
从 表 8 可 以 看 出 (最 优 值 已 加 粗 ), 在 50 次 的 函数 测试 中 ， 函数 测试 函数 
QIBAS 算法 在 对 太太 太太 所 函数 的 计算 都 收敛 到 了 最 小 ^ Rastrigin 
值 0， 标 准 差 也 为 0， 证 明了 算法 的 稳定 性 较 高 。 其 他 的 如 fo Ackley 
fr fo fos fio A. 个 函数 虽 没 有 收敛 到 最 小 值 ， 但 与 其 他 算法 相 fa Griewank 


比 其 收 全 效果 都 要 好 ,而 太 2 个 函数 与 PSO 相 比 收敛 效果 
表 8 算法 测试 结果 对 比 


Tab.8 Comparison of algorithm test results 


m GA PSO BAS QIBAS 
d 平均 值 标准 差 平均 值 标准 差 平均 值 标准 差 平均 值 标准 差 
fi 25.9571 1.3570 0.0823 0.0322 52.9531 3.9961 0 0 
f: 47.7552 3.2680 1.4682 0.5800 84.8440 3.6806 0 0 
为 8.62E+04 9.35E+03 4.47E+02 8.12E+03 9.69E+05 5.88E+06 0 0 
Ja 0.7854 0.0233 0.0455 0.8096 0.9957 3.40E-03 0 0 
f 3.499E+02 19.7682 8.65E-07 1.40E-06 7.00E+02 11.9674 2.01E+02 7.08E-07 
无 2.00E+02 0 0 0 99.80 6.4175 99.22 5.9395 
f£ 5.65E+04 3.35E+02 2.61E+03 9.5392 2.82E+03 4.30E+02 3.29E-04 3.15E-04 
fs -81.1058 1.2893 -]1.69E-02 9.13E-04 -22.7833 7.2362 -1.33E+03 7.84E+03 
为 9.91E+02 30.1911 1.5973 0.9359 -6.25E+02 58.9376 -1.99E+03 0 
à fio -4.0824 0.0231 -4.7104 3.90E-03 -20.2773 0.1413 -21.7376 1.3725 
= fi 0.1748 7.91E-03 5.62E-04 2.63E-04 1.0136 1.12E-03 0 0 
f 从 表 9 fuz ERIEUE HS, WARAPA EAP, ALORA T RA, EHRT, B 
均值 最 小 ， 由 此 得 出 QIBAS 的 性 能 强 于 其 他 三 种 算法 。 证 明了 本 文 算法 的 稳定 性 和 可 靠 性 。 
从 表 10 可 以 看 出 , 不 论 单 峰 函 数 或 者 多 峰 函数 , QIBAS 为 了 更 直观 的 展现 各 算法 在 测试 函数 上 的 收敛 情况 ， 选 
算法 和 BAS 算法 的 运行 时 间 远 小 于 PSO 算法 和 GA 算法 ， 取 结 果 与 50 次 实验 平均 值 最 为 接近 的 实验 ， 同 时 对 比 


In 


秆 间 最 高 相差 达到 100 多 倍 。QIBAS 5E BAS 两 个 算法 的 运 VSBAS SHE. SABAS 8k. BSAS 算法 以 及 BAS 算法 ， 绘 
行 时 间 相 差不多 ，QIBAS 所 需 时 间 要 长 一 些 。PSO 和 GA 各 s usc S Hl £s XE Cs Fe fos fu VU AS HU v ER cre s [8] 2E FE 
L25. EJ, QIBASh 和 BAS 运行 所 需 时 间 要 小 的 多 , 算法 K=200 进行 收敛 曲线 对 比 ,， 如 图 3~6。 E 7 Æ K 9500 HF fi. fs 


性 能 更 优 。 函数 的 收敛 曲线 。 图 8 是 K=1000 时 f. fs ERECTA ORC EL ZR 
为 了 进一步 检验 二 次 插值 对 BAS 算法 带 来 的 效果 , 分别 表 9 Friedman 检验 
进行 天 牛 须 结合 模拟 退火 (SABAS)、 变 步 长 天 牛 须 (VSBAS)、 Tab.9 Friedman test 
天 牛 群 算法 (BSAS)， 共 进行 50 次 测试， 测试 结果 如 表 11。 算法 秩 均 值 算法 秩 均 值 
K 11 的 计算 平均 值 和 标准 差 可 知 ，SABAS 算法 和 QIBAS 1.64 BAS 3.25 
VSBAS 相 比 于 标准 BAS. 算法 其 性 能 都 有 不 同 程度 的 提升 ， PSO 2.55 GA 3.49 
但 其 效果 并 不 太 明 显 。QIBAS 算法 相 比 于 SABAS 算法 除了 表 10 算法 运行 时 间 
及 大， 其 他 函数 的 优化 效果 明显 QIBAS 远 高 于 SABAS。 而 Tab.10 Algorithm running time 
相 比 于 VSBAS 算法 , QIBAS 在 fe fo 函数 的 计算 结果 稍 差 ， 算法 时间/s 大 时 间 As 时 间 /s 
CRIT. VSBAS 算法 。 大 部 分 测试 结果 优 于 BSAS QIBAS 2.860 5476 7.766 
算法 。 PSO 135.736 169.136 278.286 
从 表 12 的 检验 来 看 ， 本 文 算法 的 秩 均 值 比 SABAS、 Midi NA AM 
VSBAS 和 BSAS 等 改进 算法 都 要 小 ， 因 此 性 能 要 更 优 。 


结合 表 8 和 10 可 以 看 出 ，50 次 的 测试 结果 ， 所 得 的 标 2 :> as 
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录用 定稿 BAE, $: 基于 二 次 插值 的 天 牛 须 搜索 算法 第 38 卷 第 3 期 
表 11 算法 测试 结果 对 比 
Tab. 11 Comparison of algorithm test results 
a BAS SABAS VSBAS BSAS QIBAS 
BE E E E E E E E E FE A 
f 59.1806 3.9961 41.6044 3.9460 55.5218 3.5925 4.93E-30 0 0 0 
万 84.8440 3.6806 71.2043 3.5957 99.3674 3.9986 2.22E-15 2.6354 0 0 
f 9.69E--05 5.88E--06 9.46e-04 7.25e-04 2.47e*02 3.91e+03 2.32E-41 2.68E+02 0 0 
f 19.7657 3.40E-03 19.1583 0.2346 17.6290 1.12E-02 1.57E-15 1.0238 0 0 
Ís 7.00E+02 11.9674 45.5127 9.73E+02 1.04e+04 1.09e+03 0.9899 9.1165 2.01E+02 7.08E-07 
fs 99.80 6.4175 94.9400 6.8436 98.4400 6.0243 730 2.96E+02 99.22 5.9395 
"E 2.82E-*03 4.30E+02 1.73e+03 2.97E+02 4.09E+02 4.30E+02 3.37E-05 2.15E-02 3.29E-04 3.15E-04 
fs -22.7833 7.2362 -35.7983 6.7151 -21.4780 6.9626 -87.7179 6.8965 -1.33E+03 7.84E+03 
fs -6.25E*02 58.9376 -9.05E+02 56.9366 -88.1146 61.7779 3.9798 8.21E+02 -1.99E+03 0 
fo -20.2773 0.1413 -20.8196 0.1154 -29.5700 0.4623 -21.1639 2.3569 -21.7376 1.3725 
fa 7.5749 1.12E-03 6.2596 8.25e-04 8.0248 9.02e-04 10.2728 3.21E+02 0 0 
K 12 Friedman 检验 "m 
Tab.12 Friedman test 
算法 秩 均 值 算法 秩 均 值 
QIBAS 1.64 VSBAS 2.86 d 
BSAS 2.26 BAS 3:25 
SABAS 2.T1 -500 
从 图 3-6 可 知 ， 在 维度 K=200 时 ， 对 于 不 同 函数 BAS、 — 2 
BSAS 算法 和 SABAS 算法 的 收敛 效果 各 有 优 劣 ， 而 VSBAS -1000 
算法 的 收敛 效果 一 直 优 于 BAS 算法 , 但 与 本 文 改 进 的 算法 相 
比 , 效果 并 不 显著 。 从 图 中 可 以 看 出 QIBAS 算法 的 收敛 效果 T -ens 
远 远 高 于 其 他 四 种 算法 。 | QIBAS 
B0 r 
xi 2000 200 ECCE NEXT 4600 1800 2000 
达 代 次 数 
eo 图 5 rp fs SICHERE 
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Fig.3 Function fı convergence performance 
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图 4 函数 三 收 敛 性 能 


Fig.4 Function f, convergence performance 
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图 6 函数 fu 各 算法 收敛 性 能 对 比 
Fig.6 Function fi; convergence performance 

从 图 7 和 8 uf Al, BAS 算法 、VSBAS 算法 和 SABAS 算 
法 在 维度 分 别 为 X=500 和 天 =1000 时， 基本 没有 太 大 的 收敛 
效果 ， 甚 至 不 收敛 , 而 BSAS 算法 只 是 对 函数 的 收敛 效果 较 
好 。 这 也 证 实 了 BAS 算法 在 高 维 空间 搜索 的 局 限 性 。 但 本 文 
改进 的 算法 却 依 旧 有 较 好 的 收敛 性 , 甚至 在 X=1000 时 , 针对 
万 函数 的 收敛 效果 比 在 K =500 的 最 优 值 更 接近 实际 最 优 值 。 
为 了 更 进一步 验证 在 更 高 维 的 问题 中 ， 本 文 算法 的 可 行 
性 ， 设 计 空间 维度 X=5000,10000 ， 其 结果 如 图 9 和 10 所 示 。 
从 图 中 可 以 看 出 在 维度 为 5000，10000 时 ， 本 文 算法 依旧 有 
更 好 的 寻 优 性 能 。 
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图 7 K-500 时 函数 收敛 性 能 


Fig.7 K=500,convergence performance of functions 
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图 8 K-1000 时 函数 收敛 性 能 


Fig.8 天 =1000 , convergence performance of functions 
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图 9 K—-5000 时 各 函数 收敛 性 
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Fig.9 K =5000 , convergence performance of functions 
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综合 上 图 可 以 得 出 结论 ， 本 文 改 进 算法 (QIBAS) 在 维度 
不 断 加 高 时 ,依旧 有 较 高 的 收敛 精度 ,还 有 极 快 的 收敛 速度 。 
而 且 在 针对 不 同 函 数 时 ， 本 文 改 进 算 法 的 收敛 效果 均 较 好 ， 
也 侧面 突出 了 改进 策略 的 稳定 性 。 而 其 他 四 种 算法 却 随 着 维 
度 的 加 高 ， 针 对 不 同 函 数 来 说 ， 有 些 函 数 存 在 收敛 效果 ， 但 
对 比 低 维 度 时 效果 大 大 减弱 ， 收 敛 速度 也 大 大 降低 。 而 有 些 
则 存在 不 收敛 的 情况 。 这 也 更 进一步 验证 改进 策略 的 稳定 性 
和 有 效 性 


4 ”结束 语 


为 解决 天 牛 须 搜 索 算 法 在 高 维 空间 中 搜索 精度 低 和 易 陷 
入 局 部 最 优 的 问题 ， 本 文 提 出 一 种 基于 二 次 插值 的 天 牛 须 搜 
索 算法 .通过 二 次 插值 的 方式 提升 天 牛 跳出 局 部 最 优 的 能 力 
也 有 效 的 提升 了 天 牛 在 高 维 空间 中 的 搜索 精度 。 仿 真实 验 里 
对 11 个 测试 函数 进行 计算 ,结果 表明 ， 相 比 于 遗传 算法 、 粒 
子 群 算法 和 标准 天 牛 须 算 法 ， 改 进 算法 在 计算 精度 、 稳 定性 
和 收敛 速度 上 有 明显 的 提升 。 同 时 也 比较 了 引入 模拟 退火 的 
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Fig. 10 K =10000 , convergence performance of functions 


天 牛 须 算 法 、 变 步 长 策略 的 天 牛 须 算 法 以 及 天 牛 群 搜索 算法 ， 
统计 分 析 和 计算 结果 均 证 明 本 文 改 进 策略 的 有 效 性 和 稳定 性 
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